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We propose that the cold dark matter (CDM) is composed entirely of quark matter, arising 
from a cosmic quark-hadron transition. We show that compact gravitational objects, with 
masses around 0.5 Mq, could have evolved out of the such CDM. 

The present consensus in cosmology is that the universe is flat (D ~ 1), the baryons contribut¬ 
ing only about 10% of the total energy [i.e. VIb ~ 0.1). There is an abundance of other matter, 
the cold dark matter (CDM), which accounts for clumping on small (galactic/supergalactic 
scales), amounting to ^cdm ~ 0.35. The rest of the closure density arises from some kind of 
vacuum energy, as yet very poorly understood and termed dark energy, which is believed to be 
responsible for the accelerated expansion of the universe. The estimate of Ds comes from Big 
Bang Nucleosynthesis (BBN), whose success is one of the basic tenets of the standard cosmolog¬ 
ical model. It is thus believed that the CDM is nonbaryonic; speculations about the nature of 
CDM, all essentially beyond the standard model of particle interactions (QCD and Electroweak), 
abound and search for these exotic particles is a most active field of research. 

In recent years, there has been experimental evidence for at least one form of dark matter 
- the Massive Astrophysical Compact Halo Objects (MACHO) - detected through gravitational 
microlensing effects |3[. Based on about 13 - 17 Milky Way halo MACHOs detected in the 
direction of LMC - the Large Magellanic Cloud, a Bayesian analysis yields their mass estimate 
in the range (0.15-0.95) Mq, with the most probable value being 0.5 Mq [^,^], substantially 
higher than the fusion threshold of 0.08 Mq. (It is thus hard to understand why they do not 
ignite.) The MACHO collaboration Q suggests that the lenses are in the galactic halo. In such 
circumstances, MACHOs cannot be composed of normal baryons for reasons mentioned above. 
There have, however, been suggestions ||6|,^ that they could be primordial black holes (PBHs) 
( ~ 1 Mq ), arising from horizon scale fluctuations. This suggestion requires a fine tuning of 
the initial density perturbation and has been criticised in the literature 

Within the lore of the standard model, there occurred a phase transition from the quark- 
gluon phase to the hadronic phase during the microsecond era after the initial Big Bang, at 
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a temperature of ~ 100 MeV. The order of this phase transition is still unsettled lattice 
calculations suggest that in a pure (i.e. only gluons) SU{3) gauge theory, it is of first order. In 
the presence of dynamical quarks on the lattice, the situation is more complicated. However, 
in the early universe, the large size of the system and the long timescale could facilitate a first 
order transition. Such a transition could be modeled through a bubble nucleation scenario 


and Witten |]T| argued that the trapped false vacuum domains (TFVD) {i.e. the quark phase) 
could contain a substantial amount of baryon number. QCD-motivated studies [11,^] of baryon 
evaporation from such TFVD (called strange quark nugget or SQN hereafter) showed that if the 
SQNs contain baryon number in excess of they would be stable on cosmological time 

scales and be viable candidates for CDM, as they would be extremely non-relativistic. Note that 
the number of baryons (which would subsequently take part in BBN) within the event horizon 
at the microsecond epoch is thus the baryon number contained in SQNs could be 3-4 

times as much and SQNs of size > 10^^ could be easily accommodated. However, it has to be 
reiterated over and over again that the baryon number contained in the SQNs is in the form of 
quarks and they do not participate in BBN at all. 

We can estimate the size of the SQNs formed in the first order cosmic QCD transition in the 
manner prescribed by Kodama, Sasaki and Sato 


for details, please see Alam et al and 
Bhattacharyya et al [15|. Describing the cosmological scale factor R and the co-ordinate radius 
X in the Robertson-Walker metric through the relation 


ds'^ = —dt^ + Rfidx^ = —dt^ + R?{dX^ -|- X‘^{siv?9d(lfi + dO'^)}, 


( 1 ) 


one can solve for the evolution of the scale factor R{t) in the mixed phase of the first order 
transition. In a bubble nucleation description of the QCD transition, hadronic matter starts to 
appear as individual bubbles in the quark-gluon phase. With progressing time, they expand, 
more and more bubbles appear, coalesce and finally, when a critical fraction of the total volume 
is occupied by the hadronic phase, a continuous network of hadronic bubbles form (percolation) 
in which the quark bubbles get trapped, the TFVDs. The time at which this happens is the 
percolation time tp, whereas the time when the phase transition starts is denoted by ti. Then, 
the probability that a region of co-ordinate radius X lies entirely within the quark bubbles would 
obviously depend on the nucleation rate of the bubbles as well as the coordinate radius X{tp, U) 
of bubbles which nucleated at ti and grew till tp. For a nucleation rate I{t), this probability 
P{X,tp) is given by 


P{X, tp) = exp 


dtI{t)R^{t)[X + X{tp,ti)f 

3 Jti 


( 2 ) 


After some algebra [^], it can be shown that if all the CDM is believed to arise from SQNs, 
then their size distribution peaks, for reasonable nucleation rates, at baryon number ~ 
evidently in the stable sector. Recalling that ~ 0.1 corresponds to baryons within 

the horizon at the microsecond epoch, the total baryon number contained in SQNs to account for 
^CDM ~ 0-35 would imply 10^“® SQNs within the horizon just after the QCD phase transition. 
Because of their enormous mass, they would be nonrelativistic immediately after their formation. 
R may thus be remarked that they would be discrete macrosopic bodies (radius Rw ~ Im) 
separated by rather large distances (100 - 300m) in the background of the radiation fluid. 

Any deviation from a uniform distribution of SQNs should result in a large attractive force, 
under which they should gravitate toward one another. Given the property that they become 
more and more bound with increasing mass 0, they should tend to coalesce and grow to 
larger sizes. However, the radiation pressure acting on the moving SQNs would serve to inhibit 
such motion till such time when the gravitational force dominates over it. We can estimate the 
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Figure 1. Variation of -Fgrav/-Frad with temperature. The dot represents the point where the 
ratio assumes the value 1 . 


relative magnitude of these two forces in a straightforward manner. If the number of SQNs 
within the horizon at the time tp is 10 ® (see above), then the total number and the number 
density of SQNs at any later temperature T is given by 


Nn{T) = 10® 


100 MeVY, 

f ) ' 


nAr(T) 


Vh 


3N]^ 

47r(2t)3 


( 3 ) 


since the horizon length in the radiation dominated era is 2t. The time t and temperature T 
are related by the relation t = 0.3g* with gf* ~ 17.25 after the QCD transition [T^ ]. 

The expression for the gravitational force as a function of temperature T can be written as 
Tgrav = ^— 7 ^ where M^v is the SQN mass. (For the sake of simplicity, we assume that all 

SQNs have the same mass.) rnn{T) is the mean separation between the two nuggets, estimated 
from the density of nuggets at the temperature T. The force due to the radiation pressure on 
the SQNs arises only due to their relative motion; when they are at rest, there is no resultant 
force on them. But when the SQNs are in motion, the radiation fluid in front of the moving SQN 
gets compressed and thus exerts an additional pressure opposing the motion. In a relativistic 
framework,[| this amounts to a force F)ad given by Fj-ad = |PradCWfaii('^.RAr)/57 where prad is the 
total radiation energy density, counting all relativistic species at the temperature T, Ufaii (or 
/3c) is the velocity of the SQN and 7 the corresponding Lorentz factor. The ratio of these two 
forces, Fgrav/Frad is plotted against temperature in Fig. 1 for SQNs of initial size 10^^. Fig. 1 
readily reveals that the ratio Tgrav/Frad is very small initially. As a result, the nuggets will 

^Even though the SQNs are nonrelativistic, a relativistic treatment is necessary to handle the radiation fluid. 
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remain separated due to the radiation pressure. For temperatures lower than a critical value Tci, 
the gravitational force starts dominating, facilitating the coalescence of the SQNs under mutual 
gravity. It should come as no surprise that even for very small values of (3, the large surface area 
of the SQN is responsible for a considerable resisting force due to radiation pressure. 

We can estimate the size of the coalesced SQNs from the number of SQNs within the horizon 
at Tci- This is, of course, a lower limit, as the collapse would only start at T^i and would take 
some time during which more SQNs will enter the horizon. For SQNs of size 10^^, the total 
mass in corresponding matter within the horizon at Td turns out to be O.I 2 M 0 ; for 10^^, it 
is O.OIMq. The actual value could be much (3-10 times) higher but that can be ascertained 
only through a detailed simulation. Such a calculation is rather involved and remains a future 
project. In any case, it can be safely assumed that the coalesced SQNs will have masses above 
the fusion threshold of O.OSMq and once coalesced, the density of the resulting configuration 
would be so low that they cannot clump any further. They could thus persist till the present 
time and manifest themselves as MACHOs. For CIcdm ~ 0.35, there would be about 10^^“^“^ 
such objects within the horizon today and about 2 - 3 x 10^^ within the Milky Way halo. We 
should verify whether this abundance is consistent with the observed number of MACHO events. 

The abundance of gravitational lenses is estimated through the optical depth r which is the 
probability that a given star lies within the Einstein ring of a lens, i.e. the number density 
of the lenses times the area of the Einstein ring of each lens. The expression reads [^ ] r = 
f p(x)x(l — x)dx where Dg is the distance between the observer and the source (a star 
in the LMC in the present case) and x = D^Ds~^, being the distance between the observer 
and the lens. In particular, p is the mass-density of the MACHOs, which is of the form p = po^ 
in the spherical halo model. Assuming a halo extending all the way upto the LMC, we obtain 
r ~ 2 - 5x10“^, in excellent agreement with observation |^. 

To conclude, we have shown that in a first order cosmic quark-hadron phase transition, cold 
dark matter could arise entirely within the framework of the standard model of particle inter¬ 
actions. The observed halo MACHOs could be the natural manifestation of such CDM. 
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